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Abstract 



Stimulated by the experimental observation, made by the E760 Collabo- 
ration at Fermilab, of the ^Pi state of charmonium resonantly formed m. pp 
annihilation, we perform a calculation of the decay width for the ^Pi pp 
process. To this end, we employ a phenomenological model which adds con- 
stituent quark mass corrections to the usual massless QCD models for exclu- 
sive processes. For massless models, in fact, the process under consideration 
is forbidden by the so-called helicity selection rules, while it is allowed in our 
extended model. We find r(^Pi pp) to be in the range 1 — 10 eV. We 
also compare our results with previous, indirect estimates, based on QCD 
multipole expansion models. 

13.25.Gv, 12.38.Bx, 14.40. Gx 
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Typeset using REVT^ 



I. INTRODUCTION 



Recently the E760 Collaboration at Fermilab have reported the first possible obser- 
vation of the ^Pi bound state of charmonium (this resonance, which possibly needs further 
confirmation, has been indicated as hc{lP) in the last Review of Particle Properties 0). 
The experimental study of the ^Pi charmonium state is very interesting for testing QCD 
theoretical models. The position of its mass with respect to the center of gravity of the 
Xci^Pj) states provides a way to test different potential models for the cc bound state. Fur- 
thermore, the branching ratios for the ^Pi hadronic exclusive decays relate to the validity of 
the so-called QCD helicity selection rules (HSR) and of the QCD multipole expansion 
models (see, e.g., Ref. [Q and references therein). In particular, the ^Pi ^ pp decay is 
forbidden by the HSR's in massless QCD models for exclusive hadronic processes at high 
transfer momentum and, due to time reversal invariance, the same should hold for the 
exclusive formation process pp —>■ ^Pi, in contradiction with the results of Ref. if the 
identification of the resonance observed at Fermilab with the ^Pi state is correct. 

As a matter of fact, while massless QCD models compare reasonably well with experi- 
mental measurements in several cases, they are in disagreement with the data for a number 
of processes, particularly with those violating the HSR's. We must bear in mind, however, 
that these models and their main qualitative results, like the QCD Counting Rules or 
the HSR's, are in principle reliable at very high transfer momentum, while the available 
experimental results are all in a range of energies where their full validity can be at least 
questionable. 

Several higher twist effects can in principle be responsible for the violation of the HSR's. 
If implemented in the models, they could give, at least in some cases, a way to reconcile the- 
oretical predictions and experimental measurements. Unfortunately, a consistent inclusion 
of all higher twist effects is far from being a trivial work. In the last years, however, a few 
extensions of massless QCD models have been proposed, in order to study the possible role 
played by some of these higher twist effects. As an example, the diquark model for baryons 
is one of these (see, e.g., Ref. for a general overview on this topic). In previous work 
10-^ we have instead proposed a different, phenomenological model that extends the original 
massless QCD models by including constituent quark mass effects for the light (as compared 
to the transfer momentum scale in the process) hadrons produced in the final state. In this 
model the valence quarks of the final hadrons are given a (constituent) mass, nig = Xqirtn-, 
where Xq is the light-cone fraction of the hadron momentum carried by the quark. In the 
small-intermediate region, in fact, one may think that the valence, constituent quarks 
(the current quarks surrounded by their cloud of gg pairs and gluons) still act as a single, 
effective particle. As shown by Weinberg [0, these constituent quarks can indeed be treated 
as bare Dirac particles, with the same coupling as for current quarks in the standard model 
QCD Lagrangian. 

This model introduces spin-flip contributions in the elementary amplitudes for the hard 
scattering process, and as a consequence violations of the HSR's are in principle allowed. 
This is possibly the simplest extension of the ordinary massless models, and can be easily 
compared with them. It is assumed to be useful in the intermediate range of momentum 
transfer where, like in charmonium decays, the use of QCD models is reasonable but at 
the same time higher twist effects, of the order im?^/Q'^ {niH being a hadronic mass scale), 
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may still play an important role. As the squared momentum transfer Q"^ increases, the 
predictions of the model smoothly tend to those of the usual massless models. 

This model has been previously applied to several interesting processes , in order to 
test the role of higher twist effects at the scale of momentum transfer presently accessible 
by experiments. For example, it predicts a decay width for the HSR-violating process 
Xco VP comparable with those relative to the HSR-allowed decays Xci,c2 ^ (however, 
only a huge experimental upper limit is available for r(xco Pp), so that this prediction 
cannot be severely tested at present); it also reproduces satisfactorily the experimental 
angular distribution of the pp pair in the decay J/ip ^ pp, where the role of HSR-violating 
contributions is essential (see Ref. p[ and references therein). 

However, we must stress that other nonperturbative effects than the higher twist effects 
implemented in the model proposed here can also be relevant for charmonium decays. In 
Ref. 0, e.g., it has been shown that constituent mass corrections to the usual perturbative 
QCD models are not sufficient to explain the large experimental width for the HSR-violating 
decay rjc pp. In fact, the rjc state seems to be very peculiar, and in order to explain its, 
somehow surprisingly, large decay widths into pp or vector meson pairs, nonperturbative 
explanations, like strong gluonic components or instanton- induced contributions ||11|| , 



have been proposed. The so-called pn puzzle and the radiative decays of the J/V' |T/ 
offer other examples of processes where nonperturbative effects can play a substantial role. 

It is well possible that different higher twist or nonperturbative effects, like those quoted 
above, affect charmonium decays, each of them being perhaps more relevant for some specific 
process and less for others. A detailed comparison of different, theoretical models with all 
the available experimental results is then required to improve our present understanding of 
charmonium decays. 

Following this line of research, in this paper we apply the model of Ref. to the 

calculation of the decay width for the charmonium decay process ^Pi —>■ pp, whose interest 
has been briefiy recalled before. 

The plan of the paper is the following: in Sect. II we present the derivation of the decay 
width for the process ^Pi —>■ pp in the framework of QCD models including constituent 
quark mass effects. Since this can also be useful as a general example of how these kind of 
calculations are performed, we will give a detailed presentation of all the crucial steps of the 
calculation. In Sect. HI we present the results in a form which possibly exploit the model 
at its best, minimizing whenever possible model ambiguities. Finally, in Sect. IV we give 
some final comments and remarks. 

II. EVALUATION OF THE DECAY WIDTH FOR THE PROCESS ^P^ ^ pp 

The general expression of the differential decay width for the process ^Pi —^pp is the 
following: 

dTj^Pl^Pp) _ 1 , 4^2^1/2 V n (^P) 

d(j)d{cose) 8(27r)5 m,m^x,x' 

X A;,,;,,.m(^,0)A*_A';M'(^,0) , (l) 
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where 9 and are respectively the polar and azimuthal angles defining the outgoing direction 
of the final proton, in the ^Pi rest frame; ei is the ratio between the proton mass, m^, and 
the ^Pi mass. Mi; pi^Pi) is the spin density matrix of the ^Pi; finally, Aa,a';m(6', 0) is the 
rest frame helicity amplitude for the decay of the ^Pi bound state, with z component M of 
its total angular momentum, J = 1, into a proton and an antiproton whose helicities are 
respectively A and A'. 

It is known from first principles [|T3[ that the amplitudes Aaa';m(6', 0) have the following 
general structure: 



A 



(2) 



where the "reduced" amplitude Aw is independent of M and the angular variables and the 
d •^{6) are the usual spin rotation matrices. 

Inserting Eq. into Eq. (|l]) and integrating over dcj) we get, since /q^'^ (i0exp[i(M — 
M')<p\ = 27t6{M -M'), 



rfr(iPi ^ pp) 

d{cos9) 



K2vr) 



AA' 



A, A' 



M 



(3) 



Integrating now over d{cos6) and using the facts that Jq d{cos( 

whatever fi and fi' are and that, by parity invariance, \A p = \A^ 

we finally find: 



^ and |v4_4. 



2/(2 J+1), 
' = \A^-\\ 



TCPi^PP) 



967r4 



;i-4ef 



2U/2 



(4) 



Before proceeding with calculations, let us remind why massless QCD models predict 
a vanishing decay width for the process under consideration. First of all, amplitudes with 
the proton and antiproton having opposite helicities, Ax^-\, must vanish: in fact, the initial 
state of the process, the ^Pi, is a J^*" = 1+^ bound state; on the other hand, the outgoing 
pp pair is, for Ax^^x, in a S" = 1 state. Then, in order for the pair to have charge conjugation 
C = (—1)^"'"'^ = —1, its relative angular momentum L should be even, contradicting the 
requirement imposed by parity conservation, P = (—1)^+^, that L should be odd |]T4 . 
The crucial point is now that in massless QCD models helicity amplitudes with the proton 
and the antiproton having the same helicity, Ax^x, are also vanishing, due to the HSR's, 
which forbid spin-fiips in the elementary quark-gluon vertices. Then there is not any term 
contributing to the total decay width, Eq. (^. Notice that HSR's are valid to all orders in 
the perturbative expansion in the strong coupling constant, and as such overcoming their 
restrictions is not a matter of going to higher order in this expansion. On the contrary, 
every higher twist contribution violating the HSR's would in principle give a non vanishing 
contribution to Ax^x, of the order mn/Mi, where tjih is a typical hadronic scale. 

We now proceed considering in more details the helicity amplitudes Axx';m{0, (j)) which 
enter the calculation of the r(^Pi — ^ pp) decay width. Let us define: i) Axy'^j^i^^'P) as the 
helicity amplitude for the decay of a cc bound state with quantum numbers L, S, J, M, into 
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a pp pair with helicities A, A' respectively; ii) Mw -x^x-ik; 9, 0) as the hehcity amphtude for 
the decay of two free c, c quarks with relative momentum A; and helicities Ac, Ag into the 
same final state as before. Then, in the well-known nonrelativistic approximation, A^f,'/^ 
can be obtained by integrating the amplitude Mxx'-\^\^ over the proper cc bound state, 
nonrelativistic wave function: 



Ac Ac 

X / d^k Mxx';x^x,{k- 9, <P)Diix{P. «, 0) iJc{k) , (5) 



where the C"s are the Clebsh-Gordan coefficients required to give the right combination of 
quantum numbers; A = Ac — Ag; fc = (A;, a, (3) is the relative momentum between the c and c 
quarks and finally ipc{k) is the (momentum-space) charmonium wave function. In particular, 
for L = 1 bound states, and taking a full nonrelativistic approximation (that is, fc ^ 0), 

iljL=m = i?>V2^\R'M\^^5{k) , (6) 

where |-R'i(0)| is the absolute value of the first derivative of the bound state radial wave 
function at the origin. 

Using Eqs. (||,||,^, and defining u = cos a we get: 

A++;o(^^ = = O) = i++ = 

33/2 1 r2TT r+1 

X |Af++-++(*; 9=0=0) - Af++;__(fc; 9 = 0=0)] , (7) 



and 



A + _;l(^ = = O)=^ 

Q\3/2 1 .2tt r+1 , ,1/2 



|i?;(0)| lim - / d^ e'^ / duil 
fc^o k Jo J -I ^ 



u 



X [M++;++(fc;^^ = = O)-M++;__(fc;^^ = = O)] . (8) 

Notice that from now on we restrict our calculations to the simple kinematical config- 
uration 6^ = = 0, since the angular dependence of the helicity amplitudes Ax^x';M can be 
deduced from Eq. (||). 

In order to proceed further and evaluate the hehcity amplitudes Mxx'-XcXc{k;9,(j)) we 
resort to the previously mentioned QCD models 0. According to them and to factoriza- 
tion theorems, the amplitudes M are given, under particular conditions that are reasonably 
applicable to charmonium decays, by a convolution between a hard, perturbative contri- 
bution, describing the interaction between the cc pair and the valence quarks of the final 
baryons, and a soft, non-perturbative contribution, which accounts for the hadronization of 
the valence quarks into the final hadrons. The latter contribution is described, for each final 
hadron, by the respective distribution amplitude (DA), which is in principle independent of 
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the particular process under consideration. As is well known, the exact form of the DA's is 
essential for the quantitative results of the model. Then we have: 



Mxx'-x,x,{k,e,(f)) = J [dx][dy]ijjp^x{x)tpp,x'{y) 
X Ta,^a,2A,3,a,-,a,-2A,-3;A,a,(£, y; k, e, 0) , (9) 

where J[dz] stays for dzidz2dzs6{l—zi—Z2—zs) and we used 5 as a shorthand for {zi, Z2, z^); 
Xi{yi) represents the (light-cone) fraction of the proton (antiproton) four- momentum carried 
by the i-th quark (antiquark); the ip{z) are the baryon wave functions, containing a color, 
spin-flavor and dynamical part, the above mentioned DA. We will discuss in more details 
these wave functions in the following (see Eq. (|2^) and following comments). The helicity 
amplitude T describes the hard elementary scattering. To lowest order in the strong coupling 
constant there is only one Feynman graph contributing (apart from a set of similar terms 
obtained by permutations of the final fermionic lines, which can be accounted for by a proper 
choice of the hadron wave function). It is shown in Fig. |1|, where the kinematical notation 
is also defined. 

The hard scattering amplitude T has the following expression (we use {A} as a shorthand 
notation for the set of quark and antiquark helicities, Xg^Xq^Xq^, Xq^^Xq^Xq^): 

1 1 

r{A};A.A. = -'^^^^''''^^'glglgl^kl-mDikl-ml) 

xi?(,),.,({A})/?f,7(A„As), (10) 

where is the color factor which, once the convolution with the final hadron wave functions 
is made, takes the value = 5/(18-\/3); 

gl = M^[x,yi + {Xi - y,)hl] (11) 

is the squared four-momentum of the i-th virtual gluon, i = 1,2, 3; ki = c—gi and k2 = Qs — c 
are the four-momenta of the two virtual charm quarks (see Fig. |l|); rric is the charm quark 
mass (we assume, consistently with the adopted nonrelativistic approximation for the bound 
state, rric — Mi/2); moreover, 

{kl - ml){kl - ml) = -^{A + Bku + Ck\^) , (12) 

where k = k/Mi, and we have defined 

A = Jl \2xiyi - X., - yi + 2{xi - yife\ ; 

i=l,3 

B = (l -4ei)^^^ {(xi - t/i) [2x3^/3 - X3 - ?/3 

-F2(x3-y3)M] -(1"3)} ; 
C = -{l-4e^){x^-y,){xs-y3)■, (13) 

finally 
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(14) 



i=l 



and 



(15) 



Notice that Eq. (^ implies that the hehcity spinors for the c, c quarks are normahzed 
to one, while those for the final quarks are normalized to 2m. 
Inserting Eqs. (g),(^) into Eqs. (0),(D we find: 

33/2 ^1 

\R[{0)\cF{A7Tf [dx][dy]tlj,,4x) 



A 



++ 



—i- 



X iPpAv) 





3 p 
919293 



k=0 



^(0) ' 



(16) 



^'2 



^ \R[{0)\cF{A7r)'[dx][dy]^,AS:) 
Jo 



where we have defined (n = 0, 1) 



^2„2„2-^('?)/^'^P 

9i9293 



k=0 



^(1) 



(17) 



lim — 



u 



1— n 



1 — M 



2\'^/2 



Mf k^o kJ-i A + Bku + C^u^ 



2n 



X / ci/3 exp(m/5)F'^'''' 



and 



upup. 



(18) 



(19) 



Notice that in Eqs. (|T6D,(|T^ the strong coupling constant as has been left inside the 
convolution integral. In fact, strictly speaking, its (running) value depends on the momen- 
tum fractions Xi, Ui carried by the quarks participating to the hard scattering, as will be 
discussed in the next section. Furthermore, the color part of the baryon wave functions ipp^p 
is now included into the color factor cp. 

After some length but straightforward algebra we get 



A, 



29v^|i?;(0)|c^7r^— 9 / Mx][#]^p,+ (x)V^p,+ (y)n 



<^s{9} 



X 



B 
~A 



Mf Jo 

xm + XsVi + 2(xi - yi){x3 - y3)el 

-R(g)132 — R(q)231 



1 XiVi + {xi - ViYel A 



R{q)123 — R{q)213 + -R(g)312 " R{q)321 



X1X3 + yiys - 2{xi - yi){x3 - y3)el 
1/2 



+ (1 - 4e^) [{Xl - yi){R{q)312 - R(q)32l) " {X3 " y3){R{q)123 - %)213y 



(20) 
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A, 



1 



i2«y6|i?;(0)|c^7r^(l-4e?)i/2— [dx] [dy]i,,,+ {x)ijp,.{y) 



X 



n 



\ XiVi + {xi - ViYel A 

R(q)121 — iR 



R(q)112 — 'i'R(q)221 



+ [X2 - y2 



(g)212 



+ [x^ - y3) 



R{q)211 — 'i'R{q)122 



}■ 



(21) 



Considering now the explicit expression of the most general spin-flavor component of the 



proton wave function (see, e.g., Ref. and references therein). 



= 2A^{</.(123)Mi,A(l)lX2,-A(2)t/3,A(3) +</^(213)Mi,_A(lKA(2)t^3,A(3) 

- 2T(123)^zi,a(1)?X2,a(2)c/3,-a(3) + (1 ^ 3) + (2 ^ 3) } , 



(22) 



where is the so-called nucleon decay constant, (^{x) is the distribution amplitude, 
2T(123) = 93(132) + (^(231) (in obvious notation, (p{ijk) = ip{zi, zj, Zk)), we find that, being 
^{x,y) a generic function of x, y, 



[dx] [dy] ipp^+{x)ipp^+{y)^{x,y){R^g) 



123 



R 



(q)2i,) = ^Mle,j[dx] [dy]^{x,y) 



{<^(123)(^(213) - ^(213)^(123) + 2(^(312)T(132) 
2T(132)(^(312) - 2(^(321)T(321) + 2T(321)v3(321)} 



(23) 



r F'^ r 

J [dx] [dy] ^ljp,+ {x)^ljp,+ {y)^x, y){R(q)3i2 - R{q)32i) = ^^'^i^i J [dx] [dy] $(x, y) 

X {2(^(123)T(123) - 2T(123)(^(123) - 2(^(132)T(132) 

+ 2T(132)(^(132) - V3(321)v3(231) + ^(231)^(321)} ; (24) 



J[dx] 



^p,+ (5)V'j5,+ (y)^(5,y)(%)132 - ^(g)23l) 

X {2(/?(213)r(123) - 2r(123)y?(213) + y?(132)^(312) 
- (^(312)(^(132) - 2(^(231)T(321) + 2T(321)^(231)} 



^M',e, I [dx] [dy] $(x, y) 



(25) 



.^M', l[dx] [dy]^x,y) 



[dx] [dy] ipp^+{x)iljp-{y)^{x,y) [R{q)ii2 - iR{q)22i] 

X {(^^(123) + y?2(213) - 4X^(123) - ^^^132) + <^\312) + 4X^(132) 

+ ¥?2(321) -(/?2(231) +4r2(321)} ; (26) 

F"^ f 

[di] [dy] ^p^+{x)^p^4y)<!>{x,y) (%)i2i - ^^(9)212)) = -^^1 / [^^1 l^v] ^{x,y) 



{ - v?^(123) + v?^(213) + 4T2(123) + (^^(132) + ip^{312) - 4T^(132) 
V?2(321) + ^2(231) + 4T2(321)} ; 



(27) 



r F'^ r 

J [dx] [dy] y) {R{q)211 - iR(q)l22)) = "^^1 / [^^] M^] ^l^'^) 

X {(^^(123) - (^^(213) + 4T^(123) + ^^(132) - ^^(312) + 4T^(132) 

+ (^2(321) +(^2(^231) -4T2(321)} . (28) 

In the products of the T proton distribution amphtudes appearing in the previous 
equations, the first term is always intended to be a function of x, the second of y, while 
the number in brackets give the pertinent permutation of the zi,Z2,Z3 arguments. It is 
important to notice, as we will see in a moment, that in all the terms concerning the ^4++ 
(respectively A^_) amplitude the contribution to the convolution integral coming from the 
proton and antiproton distribution amplitudes is totally antisymmetric (symmetric) under 
the exchange of x and y. 

Using Eqs. (pOD ,(p3D-(^5D, we finally find, after some simple algebra: 



A, 



.5 ■ 2« 



33 



-IT 



^6i(i - 4ei)'/'F'Mm^, J m m n 



XiVi + {xi - ViYel 



1 



X 



X 



2xiyi -xi-yi + 2{xi - yi^el 2x^y^ - X3 - y^ + 2(^3 - y^^ej 
xi -yi X3- ys 



2xiyi -xi-yi + 2{xi - yi^ef 2x32/3 - X3 - 2/3 + 2(x3 - ys^ef 



xws + xsyi + 2(xi - 2/1) (x3 - y3)el 



ip{123)ip{213) + 2^(312)T(132) 



- 2(/3(321)r(321) + 2(^(123)T(123) - 2^(132)T(132) + (^(231)y?(321) 

- X1X3 + yiys - 2(xi - yi){x3 - y3)el 



X 



2(/?(213)T(123) + y?(132)^(312) - 2v?(231)r(321) 



+ {xi - yi) 
- {x-i - 7/3) 



2(^(123)T(123) - 2y?(132)T(132) + (^(231)(^(321) 
(^(123)(/3(213) + 2(^(312)T(132) - 2v9(321)T(321) 



(29) 



Notice that the amplitude is proportional to the ratio ei = rrip/Mi, so that, as it 
must be, it vanishes in the massless case, that is, in the (m^/Q^) ~ (m^/Q^) limit 
(remember that in our model mq = Xgirip, Xq being the light-cone fraction of the hadron 
momentum carried by the quark). 

It is also easy to verify explicitly that, as was discussed previously, the amplitude 
is vanishing, as it must be from first principles. In fact, Eqs. (^),(p6l)-(^) show that the 
argument of the integral over dx, dy is totally antisymmetric under the exchange of x and 

y- 
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III. RESULTS 



We are now equipped with all the ingredients required to perform estimates of the decay 
width for the process ^Pi pp. 

To begin with, let us recall very briefly some aspects of QCD models (both in the massless 
and massive cases) which are essential in determining their quantitative predictions. As 
previously mentioned, a very important component of the calculations are the DA's. The 
so-called asymptotic DA's, as deducible from perturbative QCD at very large values of the 
transfer momentum, are not guaranteed to be reliable in the range of energies experimentally 
accessible today and presumably also in the near future, due to the their slow, logarithmic 
evolution with Q^. Non-perturbative methods, like Lattice QCD or QCD sum rules, give 
information on the first moments of the DA's that lead to a number of models, widely used 
for practical calculations. We will consider several of these model DA's, in order to get 
indications on the spread of our results due to this ingredient, without giving prominence 
to one model over the others. 

The value of the first derivative of the nonrelativistic, radial bound state wave function 
is another parameter to be fixed. One may use potential models to estimate this quantity 
or, alternatively, fix its value phenomenologically. This can be done by comparing the 
theoretical prediction with the experimental measurement for some other process which is 
better understood theoretically and more accessible experimentally. 

Sudakov factors can also change sizably the quantitative predictions. The inclusion of 
these factors has been advocated in order to improve the consistency of the theoretical 
models for exclusive processes and answer some criticisms connected with the end-point 



contributions (see Ref. [16] and references therein). 

Furthermore, we do not know whether and how the DA's are modified due to higher 
twist, mass correction effects. 

In order to improve the significance of our results and, where possible, to minimize their 
dependence on all the abovementioned uncertainties, we choose to evaluate, rather than 
the absolute value of the T(^Pi ^ pp) decay width, its ratio with the decay width for the 
analogous process, Xc2 —^pp- In other words, we write 



r(ipi ^pp) 



r(ipi -^pp) 
r(xc2 ^pp) 



r(xc2^PP), (30) 



where we use the theoretical prediction for T{xc2 ^ pp) in the denominator and the corre- 
sponding experimental measurement in the numerator of the r.h.s. of Eq. (|30D . 

This way, we can neglect the dependence of the results on F/v and |-R']^(0)|. Notice that F/v 
and |-R'i(0)| appear in the expression of T{}Pi — > pp) to the 4th and 2nd power respectively, 
so that even small variations in their values might give relatively sizable changes in our 
observable. Possible Sudakov factor effects might also be mitigated in the ratio. We are 
then confident that our numerical estimates for the T{}Pi —>■ pp) decay width, given in the 
form of Eq. (pOD, have the higher precision attainable with present models. 

The decay width T{xc2 pp) has been evaluated previously in massless QCD models 
(see Ref. [jT5[ and references therein). We have done the calculation following the same 



procedure previously described for the T(}Pi ^ pp) case. Our result is in agreement with 
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that of Ref. |T5|. Notice that this calculation accounts only for the leading term in a 



{rrip/M) power expansion (in that case the {irip/M)^ term, since the process is not forbidden 
by HSR's). For consistency, we also keep only the leading, (m^/M), term in the calculation 
of T{}Pi pp). It is important to notice that, this way, we avoid also possible problems 
connected with modifications of the DA's when mass corrections are taken into account, 
since these modifications do not contribute at leading order. 
Using Eq. (W), and Eq. at leading order in ei we have: 



tCp.^pp) 



xF^|i?;(o)p 



b^n^el (l - Ael 
1 



1/2 



Ml 



12 



It 



(31) 



where Ii is the integral appearing in Eq. (^), taken in the limit ei — 0. Notice that 
we have retained the mass contribution in the kinematical factor coming from phase space 

1 /2 

integration, (1 — 4e^) 



The corresponding expression for the Xc2 pp decay is [p^ , p!^] 

2" 



r(xc2 ->pp) 



36-5 

X ^^l^'l 



TT 



4e; 



1/2 



(32) 



where M2 is the mass of the Xc2 and €2 = mp/M2. For completeness, we report the expression 



of the integral I2 (see also Ref. [p!5|] ) 



I2 = j[dx][dy]as[x^yiM']as[x^y'iM']as[{l-xi){l-yi)M'] 

xi +yi 

xiyix^yz{l - - yi) [xi{l - yi) + yi{l - xi)] 
X {^^(231) + ^2(132) + ^^(123) + (/?2(321) + 4X^(123) + 4r2(132)} . 



(33) 



Inserting Eq.s (|3T| ) ,(l3^ ) into Eq. (|30D we find finally: 



T{^Pi^pp) 



24 ■ 3 



1 -4e^ 

1^2, 



2\ 1/2 



M2 

Ml 



12 



pp) . 



(34) 



In Table | we present the results obtained using Eq. (|^) and the available models for 
the DA's. Besides the nonrelativistic (NR) and asymptotic (AS) DA's 0, we have used 
model DA's inspired by QCD sum rules methods: the Chernyak, Zhitnitsky (CZ) model 



and its improved (COZ) version, due to Chernyak, Ogloblin, and Zhitnitsky (see Ref. |T5 



and references therein); the King, Sachrajda (KS) model |]T8|; the Gari, Stefanis (GS) model 
||19|| ; and, finally, the Heterotic (HET) model proposed by Stefanis and Bergmann |2^. A 
collection of the explicit expressions of these DA's can be found in Ref. P]. 

The expression used for the strong, running coupling constant inside the convolution 
integrals /i and I2 is the following: 
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as{Q') = T > (35) 

(lln, - 2n/) log [(g2 + 4m2)/A2] 

where tt-c = 3 is the number of colors; rif is the number of active flavors (in this context, 
Hf = 4); TUg ~ 0.5 GeV is an effective, dynamical gluon mass regularizing the a.^ behavior 
at slow Q'^ ||21[; we assume also A = 0.2 GeV/c. 

The evolution of the proton DA's, according to perturbative QCD, has been also 
taken into account. 



Finally, we have used r^^P(xc2 ^ pp) = (206 ± 22) eV 0j2| . 

It is evident from Table | that, in spite of our efforts in order to limit the dependence of 
the numerical results on the particular shape of the model DA, our estimate for T{^Pi pp), 
even excluding the NR and AS distribution amplitudes, still cover a range from 1 to 10 eV. 
This spread of the results is in part due to the different behavior of the integrals Ji and 
I2 when changing the DA's. Relative changes are more sensitive for /i than for I2. This 
can be easily understood from a qualitative point of view. In fact, the integral Ji is exactly 
vanishing for DA's if{z) that, like the nonrelativistic and the asymptotic ones, are completely 
symmetric under permutations of zi,Z2,zs. Moreover, due to the factors {xi — Ui) which are 
present in each contribution to Ji, Ji itself is very depressed when ~ as is the case for 
DA's that are strongly peaked at Xj ~ ?/j ~ 1/3. Then, the magnitude of /i should give a 
measure of how much the corresponding DA is far from the simple, intuitive equipartition 
of the hadron momentum among its valence constituents. In principle this could give an 
interesting way to discriminate at least the gross features of the proton DA's. In practice, 
however, at this stage of the calculations, these effects are probably masked by other, less 
controlled, sources of uncertainty. From this point of view it should surely be of relevance 
to perform a full implementation of mass corrections. We have not calculated them for the 
Xc2 pp process, while for the ^Pi case, direct use of the full expression of Eq. (p9D, 
neglecting possible modifications in the shape of the DA, gives a sizable reduction in the 
total width. The worst case then would be if a full account of mass corrections in the Xc2 case 
should give an enhancement of the corresponding decay width, depressing our prediction for 
the relative T{}Pi pp) decay width, Eq. (^). We must not forget, however, that the 
results presented in Table |I[ which are at leading order in the mass corrections, are on the 
contrary free from the ambiguity connected with possible changes in the shape of the DA's 
due to higher order mass corrections, since they do not contribute to leading order. 

Let us now compare our predictions for T{^Pi pp) with previous estimates. Assuming 
the resonance observed by the E760 Collaboration to be the ^Pi state of charmonium, we have 
the following measurement for the product of the branching ratios B(}Pi pp)B{^Pi 
J/^Pn^) i: 

(2.3 ± 0.6) X 10"^ > BCPi pP)BQPi ^ J/^l^n^) 

> (1.7±0.4) X 10"^ (36) 

The range of values in Eq. ( p6D corresponds to the plausible range of values for the ^Pi 
total decay width, rr(^Pi) 0: 

500 keV < Vt{^Pi) < 1000 keV . (37) 
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By using Eq. (|37D and the biggest of our estimates from Table I, we find correspondingly 

2.1 X 10~^ > BCPi VP) > 1-0 X 10"^ . (38) 

Unfortunately we are not able to give a reliable prediction for the decay width of the 
process ^Pi — > J/ipir^ in the framework of our model. Then, we cannot make a full com- 
parison of our results with the experimental ones of Ref. [|I|. A complementary situation 
is found in the case of the QCD multipole expansion models. In this framework, Kuang et 
al. [Q have predicted T(^Pi J/ipn^) = 0.3 (a^^/a^) keV, where and a^^- are effective 
electric and magnetic multipole-expansion coupling constants. There is some uncertainty 
on the value of the ratio a^^/a^, but theoretical considerations suggest as reasonable range 
1 < a^^/a^ < 3 [Q. We thus take as indicative estimate a^j^/a^ ~ 2, which leads to the pre- 
diction r(^Pi J/ipn'^) ~ 0.6 keV. Comparing this result with Eqs. (|36D , (|37|) , the following 
estimate for B{}Pi ^ vP) follows: 

1.9 X 10"^ < Bi^Pi VP) < 2.8 X 10"^ . (39) 

The two predictions of Eqs. (^),(^) differ approximately of one order of magnitude; they 
are closer the smaller is the total decay width for the ^Pi state. However, both theoretical 
calculations are at present order of magnitude estimates, and before drawing any definite 
conclusion better calculations are to be expected, together with improved experimental 
results. 

Let us finally remark that our prediction for the branching ratio B{}Pi —* vP)j Eq. (PHD, 
compares reasonably well with the experimental branching ratios for the processes Xc2,i 
VP, B{xc2 VP) = (10.0 ± 1.0) X 10-^ B{xci VP) = (8-6 ± 1.2) x 10"^ In fact, in 



QCD models the ^Pi — > vP decay should be suppressed with respect to the Xc2,i VP ones 
by a factor ~ iTi'jj/Q'^, being mn ~ 1 GeV a typical mass scale for higher twist effects. For 
Q2 ^ 2 this should correspond to a suppression factor of about 10 — 15, which should in 
turn put the expected value for B(^Pi vP) the range covered by our prediction. 



IV. CONCLUSIONS 

In this paper we have estimated the width for the decay process of the ^Pi charmonium 
state into a proton-antiproton pair, using QCD models and including constituent quark mass 
corrections. In fact, massless QCD models state that this decay width is vanishing, due to 
the helicity selection rules. This is in possible contradiction with the recent experimental 
findings of the E760 Collaboration. The inclusion of higher twist effects, like constituent 
quark mass effects, which can still play a non negligible role at the involved energies, can 
in principle improve considerably the situation. However, a full, consistent treatment of 
these effects is out of our present possibilities. We have then used a phenomenological 
model which extends the massless QCD ones. In order to minimize the dependence of 
the quantitative results of the calculation on details of the model, we have expressed the 
^Pi vP decay width in connection with that for the Xc2 VP decay process: T(^Pi — *■ 
VP) = [r(iPi ^ VP)/^{Xc2 VP)Y^ X r'=^P(xc2 ^ VP)- That is, we evaluate the first 
term in this expression in our model, at leading order in the parameter e = nip/M, while 
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taking the available experimental measurement for the second term. This way, not only 
some phenomenological parameters of the model (the nucleon decay constant F/v, the value 
of the first derivative of the radial, L = 1, nonrelativistic charmonium bound-state wave 
function at the origin) cancel out from the result, but some other effects should be at 
least minimized. For example, the possible effect of Sudakov factors should be reduced. 
Furthermore, a possible dependence of the DA's on mass effects is also irrelevant as long as 
only the leading term in the e = rUp/M power expansion is kept into account. A consistent 
inclusion of non-leading terms in the e power expansion would in fact require an estimate 
of mass effects on the T{xc2 —>■ pp) decay width and on the DA's. As shown in table |, 
we estimate the Ti^Pi pp) decay width to be of the order 1 — 10 eV, depending on 
the DA considered. The corresponding branching ratio depends of course on the ^Pi total 
decay width, TTi^Pi), which is not known at present. Using, as suggested by the E760 
Collaboration, the range 500 keV < Tt(^Pi) < 1000 keV, and the biggest of our estimates of 
Table I, we find 2.1 x 10~^ > B(^Pi pp) > 1.0 x 10~^. This result is sizably lower (about 
one order of magnitude) than the corresponding branching ratios for the Xci,c2 states, as 
measured by the E760 Collaboration. This is reasonable, since in QCD models the decay 
^Pi pp should be suppressed by a factor m|^(~ lGeV^)/M^ ~ 1/12 with respect to the 
decays Xci,2 pp- The above result is also sizably lower than that obtained by comparing 
the experimental results of Ref. |]l| with theoretical estimates of the decay width for the 
process ^Pi J/ip + tt^, obtained in the framework of QCD multipole expansion models. 
Both our calculation and that based on QCD multiple expansion models can probably be 
affected by "theoretical" errors which can modify their quantitative predictions by a sizable 
factor. The two estimates, at least indicatively, can be taken as representative of two classes 
of different models: the first extending the usual perturbative QCD models by the inclusion 
of higher twist effects, like mass corrections; the latter including effectively some different, 
nonperturbative effects, like those present, e.g., in the rjc — * pp decay process. As such, their 
results should at least suggest a definite (albeit still large) range of possible values for the 
r(^Pi —>■ pp) decay width, ranging approximately from some eV to some hundreds of eV. 

We conclude observing once more that further, more precise, experimental measurements 
will certainly help in clarifying our understanding of the ^Pi —>■ pp decay process and will 
be very useful to improve in general our theoretical models for the calculation of exclusive 
decay processes in the charmonium family. 
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TABLES 



TABLE L Predictions, as from Eq. (^), for the decay width relative to the process ^Pi pp, 
considering several available models for the nucleon distribution amplitude (see text for more 
details). 



DA Tj^Pi^pp) (eV) 

l^R 00 

AS 0.0 

CZ 10.5 

COZ 2.7 

KS 1.2 

GS 1.3 

HET 4.7 



16 



FIGURES 



FIG. 1. The Feynman diagram which, to lowest order in a^, describes the elementary process 
QQ ~^ Q1Q2Q3Q1Q2Q3, for a quarkonium state with charge conjugation C = —1. In the QQ cen- 
ter-of-mass frame, = (E, k/2) and = {E, —k/2), where k is the relative momentum between 
the c and c quarks; qi = Xip and qi = yip (i = 1,2,3), with = (E,q), = {E,—q) and 
q = (gsin^ cos 0, gsin^sin^, g cos a,b,c, i,j,l,l',mi^2,3, ^1,2,3 are color indices; the A's label 
helicities. 



17 




'91 

nil 



92 

m2 



*93 



Ql^Ki <l2,K]2 (l2>i\: 



'93 



Fig. 1 



